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ABSTRACT 

Context. Apart from its contribution to cosmology, the WMAP data brings new information on the Galactic interstellar 
medium. In particular the polarization data provide constraints on the Galactic magnetic field and the synchrotron 
emission, while the intensity data can be used to study the anomalous microwave emission. 

Aims. The main goals of this study is to use the information from both WMAP intensity and polarization data to do 
a separation of the Galactic components, with a focus on the synchrotron and anomalous emissions. 
Methods. Our analysis is made at 23 GHz where the signal-to-noise ratio is the highest and the estimate of the CMB 
map is less critical. Our estimate of the synchrotron intensity is based on an extrapolation of the Haslam 408 MHz 
data with a spatially varying spectral index constrained by the WMAP 23 GHz polarization data and a bi-symmetrical 
spiral model of the galactic magnetic field with a turbulent part following a -5/3 power law spectrum. 
Results. The 23 GHz polarization data are found to be compatible with a magnetic field with a pitch angle p = —8.5° 
and an amplitude of the turbulent part of the magnetic field 0.57 times the local value of the field, in agreement with 
what is found using rotation measures of pulsars and polarized extinction by dust. The synchrotron spectral index 
between 408 MHz and 23 GHz obtained from polarization data and our model of the magnetic field has a mean value 
of f3s — —3.00 with a limited spatial variation with a standard deviation of 0.06. When thermal dust, free-free and 
synchrotron are removed from the WMAP intensity data, the residual anomalous emission is highly correlated with 
thermal dust emission with a spectrum in agreement with spinning dust models. 

Conclusions. Considering a classical model of the large scale Galactic magnetic field, we show that the polarization 
data of WMAP are in favor of a soft synchrotron intensity highly correlated with the 408 MHz data. Furthermore 
the combination of the WMAP polarization and intensity data brings strong evidence for the presence of unpolarized 
spinning dust emission in the 20-60 GHz range. In preparation for the Planck mission this joint analysis of polarization 
and intensity data opens new perspective on the study of the Galactic interstellar medium and on the component 
separation exercise. 

Key words. Polarization: Radiation mechanisms: non-thermal: Turbulence: Magnetic fields: Dust 



1. Introduction 

One of the greatest challenges of observing the Cosmic 

Microwave Background (CMB) in the 20-200 GHz range 
resides in the separation between the CMB and foreground 
emission. Tliis task is facilitated by the fact that the inten- 
sity of the emission from the Galactic interstellar medium 
reaches a minimum in this range. On the other hand, even 
if the Galactic emission is weak, it is still stronger than 
the CMB over a significant fraction of the sky. In addition 
the identification of the cosmological signal is complicated 
by the fact that several interstellar emissions arc superim- 
posed in this frequency range: free-free, synchrotron, ther- 
mal dust. There are also strong evidence of an excess of 
emission in this range, discovered by Kogut et al. (1996b, a) 
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in the COBE-DMR data. This emission, is remarkably well 

correlated with thermal dust emission at 100 /im with a 
spectrum which can be confused with that of free-free or 
synchrotron emissions. 

Leitch et al. (1997) showed that the correlation of this 

excess with 100 /im extends to sub-degree scales in the 
North Celestial Loop (NCL) region. The lack of correla- 
tion of the anomalous emission with Ha measurements, also 
highlighted by de Oliveira-Costa et al. (2002) using Tenerife 
and COBE-DMR data, led Leitch et al. (1997) to conclude 
that free-free emission can only account for the data if 
it comes from gas with Tg > 10*^ K. Draine & Lazarian 
(1998a) ruled out this interpretation based on energetic 
arguments and suggested rotational or spinning dust as a 
plausible emission mechanism (Draine & Lazarian, 1998b). 
This interpretation received strong support from dedicated 
ground observations (e.g Watson et al., 2005; Fernandez- 
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Fig. 1. Free- Free emission at 23 GHz (in logio(mK)) estimated 
using Ha emission corrected for extinction (following Dickinson 
et al. (2003)) and the WMAP MEM decomposition. For E{B - 
V) > 2, the free-free emission is the one obtained from the 
WMAP MEM decomposition. For E{B -V)<2, the free-free 
emission is the one estimated using Ha emission, expect when 
the WMAP MEM free-free is lower. 



Cerezo et al., 2006) and from the analysis of the WMAP 
data (Lagache, 2003; Finkbeiner, 2004; Davies et al., 2006). 

The WMAP team proposed an alternative explanation. 
Bennett et al. (2003a) and Hinshaw et al. (2007) argued 
that high-energy electrons, responsible for the synchrotron 
emission at WMAP frequencies, are produced in star form- 
ing regions which are also bright in the infrared, suggesting 
a natural correlation between synchrotron and dust emis- 
sion. As the synchrotron emission is less correlated with 
dust emission at lower frequencies (e.g. 408 MHz), the dust 
correlated emission observed in the WMAP data is ex- 
plained as a change in the synchrotron spectral index which 
can only be attributed to a local change of the energy distri- 
bution of cosmic rays. Spatial variations of the correlation 
between the magnetic field intensity and the gas density 
would not produce a frequency dependent behavior (see 
§ 3.4). This interpretation has one important drawback: the 
correlation between the 23 GHz WMAP emission and the 
100 /im emission extends down to low column density cir- 
rus clouds, where there is no star forming activities and no 
local production of cosmic rays (see Davies et al. (2006)). 

In this paper we want to put some constraints on the 
presence and properties of the anomalous emission in the 
WMAP frequency range, while making minimum assump- 
tions about the emission mechanisms. The difficulty of this 
task comes from our lack of knowledge on the other compo- 
nents at these frequencies, mostly synchrotron. Anomalous 
emission has been measured to be very weakly polarized 
by Battistelli et al. (2006). This can be understood if the 
carriers of the spinning dust emission are small dust grains 
which are known not to produce any polarization from ex- 
tinction measurements (e.g. Martin, 2007). Free-free is also 
known to produce no polarization^. We therefore follow 
Finkbeiner (2004) and use polarization data to estimate the 
synchrotron emission. Based on this assumption and on a 
model of the Galactic magnetic field, we use the WMAP 
23 GHz polarization data to put realistic constraints on 
the synchrotron intensity and isolate the anomalous emis- 
sion component. 

The paper is organized as follow: the WMAP data we 
used are summarized in § 2. Then § 3 describes the mod- 
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els we used for the Galactic components, § 4 describes 
our model for the synchrotron polarized data and for the 
Galactic magnetic field, and § 5 presents the anomalous 
emission extracted from the WMAP data. The results are 
discussed and summarized in § 6. 

2. WMAP data 

In this analysis we use the three-year WMAP products 
available in the Healpix pixelisation scheme (Gorski et al., 
2005) on the LAMBDA website^. We use the temperature 
and polarization (I, Q, U) coadded maps per frequency in 
the K (23 GHz), Ka (33 GHz), Q (41 GHz), V (61 GHz) 
and W (94 GHz) bands. The original angular resolution of 
each map ranges from 0.23° at 94 GHz to 0.93° at 23 GHz. 
The maps were all smoothed to a common resolution of 
1°. We also use the Galactic Kp2 and point sources masks 
provided by the WMAP team. 

To remove the Cosmic Microwave Background (CMB) 
from the data we use the estimate based on the Internal 
Linear Combination (ILC) method proposed by the WMAP 
team and applied on the three years data (Bennett et al., 
2003a). The resolution of the ILC map is also 1°. As shown 
by Davies et al. (2006); Dobler & Finkbeiner (2007) the 
estimate of the CMB can have a significant impact on the 
determination of the foreground spectra. In our case, it is a 
sub-dominant effect as our analysis is conducted at 23 GHz 
where the Galaxy-to-CMB ratio is the highest. 

By design WMAP is unable to measure the zero level 
at each frequency. This has no consequence for the study of 
the CMB fluctuations but it has an importance for the de- 
termination of the foreground properties. In practice there 
is a degeneracy between the zero level and the absolute 
foreground component amplitudes one can deduced from 
the data, unless one relies on a correlation analysis like 
Davies et al. (2006) who measured a mean spectral index 
of Ps ~ -3 between 408 MHz and 23 GHz at high Galactic 
latitude. In their analysis Eriksen et al. (2007) fitted jointly 
the CMB, foregrounds and zero level. In order to estimate 
the zero level they put a constrain on the average value of 
(3s based on the results of Davies et al. (2006). In our anal- 
ysis we rely on the zero level determined by Eriksen et al. 
(2007). In practice this zero level correction is only relevant 
at 23 GHz. 

3. Galactic components 

3.1. WMAP data model 

In the following we consider that the WMAP intensity 
(corrected for CMB, zero level and dipole, and smoothed at 
1°) at frequency v is the sum of free- free (F^), synchrotron 
{S,y), thermal dust (Di,) and anomalous emissions (Ai,)^: 

I, ^ + + + [A,]. (1) 

In polarization, the model of the WMAP data involves 
less components as the free-free and anomalous emissions 
are assumed to be unpolarized. The CMB polarization level 
is so weak (Page et al., 2007) that it is not included in the 
data model. The polarized emission model is thus simply: 

P, = SP + DP. (2) 
^ http://lambda.gsfc.nasa.org 

^ in one model (model #1) we assume no anomalous emission 
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Fig. 2. Synchrotron spectral index (/3s) between 408 MHz and 
23 GHz. All maps were smoothed at 5°. Top: Model 1 - this 
is making the assumption of no anomalous emission. The map 
is similar to the one of Bennett et al. (2003b) (their Fig. 4). 
Middle: Model 3 - assumes a spinning dust component at 23 
GHz which scales with E(B-V). Bottom: Model 4 - [3s obtained 
from polarization data for p = —8.5°, xo = 8°, vq — 11°, and 
Cturb = 1.7 ^G. 



3.2. Free- Free 

Free-free emission is a rather difficult component of the in- 
terstellar emission to identify as it does not dominate at any 
frequency. Ha, corrected for dust extinction, is often used 
as a proxy for free-free as they both depend linearly on the 
emission measure, EM= / nl dl (see Dickinson et al. (2003) 
and Finkbeiner (2003)). On the other hand, several system- 
atic errors can affect the estimate of free-free emission from 
Ha corrected for extinction (see Bennett et al. (2003a) and 
Finkbeiner (2004)). Finkbeiner (2004) estimated that an 
error of a factor of 2 is likely in some parts of the sky solely 
due to an imperfect extinction correction. Electronic tem- 
perature variations are also expected, with lower values in 
the Galactic plane due to a higher metallicity and therefore 
more efficient cooling. Estimate of spatial variations of 




Fig. 3. Histogram of (3^ between 408 MHz and 23 GHz 
smoothed at 5°. solid line: model 1, Ps — —2.83 ± 0.11 dot- 
ted line: model 3, /3s = —2.96 ± 0.07 dashed line: model 4, 
I3s = -3.00 ±0.06. 
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Fig. 4. Residual 23 GHz emission when free-free and syn- 
chrotron with a constant /3s = —3 are subtracted, as a function 
of E(B-V). The sohd line is 0.65 E(B-V) + 0.03. 



are not available and a constant value on the sky is usually 
assumed. 

It is important to point out that the free-free spectral 
index is a slowly varying function of frequency and electron 
temperature (Bennett et al., 1992; Dickinson et al., 2003; 
Bennett et al., 2003b): 



1 



10.48 + 1.51n(Te/8000AO - ln{iyGHz) 



(3) 



At 23 GHz, the free- free spectral index varies only from 
2.13 to 2.16 for 4000 < < 10000 K and it is not affected 
by extinction. Based on the analysis by Davies et al. (2006) 
we also note that, at this frequency, the spectral index of 
free-free is significantly different than that of synchrotron 
(/3 ~ —3) and anomalous emission (/? ^ —2.85). This is 
a favorable configuration for a maximum entropy method 
(MEM) analysis that can easily isolate a component over 
the whole sky that has a fixed spectral index using five 
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frequency bands, especially if the other components have 
significantly different spectral indices. 

Based on this physical property of the free-free emission 
Bennett ct al. (2003b); Hinshaw ct al. (2007) estimated 
the frcc-frcc cniission in the WMAP data using MEM and 
considering a constant spectral index of 2.14 over the whole 
sky. In the following, we rely on the result of the MEM 
decomposition of Bennett et al. (2003b) for the estimate of 
the free-free emission in regions where E(B — V) > 2 (or 
■^v > 6). The same assumption guided the work of Sun 
et al. (2007). In more diffuse regions we use the free- free 
estimate of Dickinson et al. (2003) based on a Ha all-sky 
template at 1° resolution, unless the WMAP MEM free- 
free is lower than the Ha estimate. The free-free template 
we use at 23 GHz is shown in Fig. 1. 

3.3. Thermal dust 

To estimate thermal dust emission in the WMAP data we 
used model #7 of Finkbeiner et al. (1999), multiplied by 1.2 
to account for the emission in the 94 GHz band (Hinshaw 
et al., 2007). The correlation of the predicted thermal dust 
emission with (194 — F94) at 94 GHz and inside the Kp2 
mask is 0.96. We note that this component contributes only 
to 1% of the emission at 23 GHz and that it does not 
affect our analysis of the synchrotron and anomalous emis- 
sion at this frequency. This is true also in polarization as the 
polarization fraction of thermal dust is at the level of only 
a few percent (Benoit et al., 2004; Ponthieu et al., 2005; 
Page et al., 2007). Therefore at 23 GHz, where most of our 
analysis is done, we considered the polarized intensity to 
be entirely due to synchrotron emission: P23 = 823- 

3.4. Synchrotron 

The synchrotron intensity depends on the cosmic ray (elec- 
tron) density (rig) and the strength of the magnetic field 
perpendicular to the line of sight {B±). For a cosmic ray 
distribution following a power law, N{E) cx E^'^, the syn- 
chrotron intensity at frequency v is: 



{l+s)/2 



dz 



(4) 



where the integral is over the line of sight z and B± = 

'bI^I with the plane x - y corresponding to the plane 

of the sky. The emissivity term es{i') is given by a power 
law: 



where 

/3« = -(s + 3)/2. 



(5) 



(6) 



It follows that the synchrotron intensity spectrum is a 
power law 



(7) 



For a typical value of cosmic ray spectrum s = 3, we obtain 
/3« = -3. ^ ^ 

The most reliable estimate of synchrotron emission to 
date is the 408 MHz all-sky map of Haslam et al. (1982). In 



the following we estimate the synchrotron emission at any 
WMAP frequency using an extrapolation of the Haslam 
408 MHz map (5*408) corrected for free-free emission as pro- 
vided by the WMAP team: 



'408 



0.408 



(8) 



The extrapolation of the 408 MHz emission at WMAP 
frequencies relies on an estimate of fig and of its variations 
on the sky. All-sky estimates of the synchrotron spectral in- 
dex between 408 MHz and 2326 MHz have been attempted 
by Giardino et al. (2002); Platania et al. (2003). They ob- 
tained compatible values of fig = —2.7 ± 0.1. On the other 
hand the spectral index varies with frequency (Banday & 
Wolfendale, 1991) due to electron aging effects. This effect 
is expected to be important in the WMAP range which 
makes uncertain the extrapolation to 23 GHz of the spec- 
tral index defined at lower frequencies. 

Prom the analysis of the WMAP data Finkbeiner (2004) 
found that a constant (ig = —3.05 between 408 MHz and 
23 GHz corrects the most obvious synchrotron features seen 
both at 408 MHz and 23 GHz (i.e. Loop 1). This is in accor- 
dance with the work of Davies et al. (2006) who obtained 
13s = -3.01 ± 0.04 outside the Kp2 mask. 

In the following we describe four different models of the 
synchrotron spectral index between 408 MHz and 23 GHz. 
In all models we considered that thermal dust emission is 
negligeable at 23 GHz {D2?, ^ 0) and that the free-free {F^) 
is given by our model (see § 3.2). 

3.4.1. Model 1 

Following Bennett et al. (2003a), this model assumes that 
the 23 GHz emission is only synchrotron and free-free, giv- 
ing an upper limit on the synchrotron spectral index: 



Pg = 0.248 log 



I23 — F', 



23 



'408 



(9) 



where the multiplicative factor 0.248 is simply 
l/log(23/0.408). The map of (3g obtained is shown in 
the top panel of Fig. 2. Using this simple method, which 
does not take into account anomalous emission, we obtain 
(3g = -2.84 ± 0.12. The histogram of (3g is shown in Fig. 3 
(solid line). Taking into account anomalous emission will 
lower the contribution of synchrotron at 23 GHz which 
will increase the value of /Jg- This (3s is similar to the one 
obtained by Bennett et al. (2003a). The differences are 
mainly due to the zero level correction we applied and to 
our slightly different free-free estimate. 

3.4.2. Model 2 

In this model we assume the presence of anomalous emis- 
sion making no hypothesis on its spatial distribution and 
we set a constrain on j3s considering it constant with a value 
of (3s = —3 in accordance with Davies et al. (2006): 



.^23 = Si 



'408 



23 
0.408 



+ i^23 + A 



23- 



(10) 



When this synchrotron estimate and our free-free model are 

subtracted, the residual (A23) is shown to be well correlated 
to the estimate of E(B-V) of Schlegel et al. (1998) smoothed 
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at 1 degree. In diffuse regions (E(B-V)< 0.5) we obtained 
^23 ~ (0.65 ± 0.01) E(B-V) +(0.03 ± 0.01) (see Fig. 4). 

Even though the estimate of E(B-V) of Schlegel et al. 
(1998) has been shown to overestimate systematically the 
extinction in regions where colder dust is present (see 
Cambresy et al., 2001, 2005), it remains a relatively ac- 
curate tracer of dust column density in diffuse regions, typ- 
ically where E(B-V)< 0.5. 



3.4.3. Model 3 

We can improve the estimate of Pg and estimate its spatial 
variation by assuming a template for the anomalous emis- 
sion component. Here we make the assumption that the 
anomalous emission is well correlated with dust column 
density (or E(B-V)). We added a component that scales 
with E(B-V): /23 = + ^23 + aE{B - V). Here the syn- 
chrotron spectral index is: 



and 



= 0.248 log 



hs-F2i-aE{B-V) 



'408 



(11) 



In this simple model we use a = 0.65 mK/mag which is 
the slope observed in diffuse regions in model 2 (see Fig. 4) . 
Compare to model 1, this model (see Fig. 2, middle panel 
and Fig. 3 (dotted line)) produces a steeper fig on average 
(/3s = —2.96 ± 0.07) but compatible with model 2 (con- 
stant Ps = —3) which is expected as we used the value of a 
deduced from model 2. Most importantly the spatial varia- 
tions of j3s is reduced by a factor of 2 compared to model 1. 
This has to be appreciated considering that some spatial 
variations of /3s could be due to spatial variations of the 
anomalous emission that we considered to be strictly pro- 
portional to E{B — V). Nevertheless the fact that taking 
into account an extra component reduces significantly the 
variations of /3s is a strong indication of a better modeling 
of the data. This is compatible with gamma ray observa- 
tions from whic;h Strong et al. (2000) do not predict strong 
spatial variations of the synchrotron spectral index, mostly 
due to the large interaction length of cosmic ray electrons. 

3.4.4. Model 4 

In order to separate the synchrotron and anomalous emis- 
sions at 23 GHz we would like to have an estimate of /3s 

independent of intensity measurements at 23 GHz and of 
any constrains on the spatial variation of the anomalous 
emission. In model 4 we use the fact that the polarized in- 
tensity at 23 GHz is dominated by synchrotron. Based on a 
model of the Galac;tic; magnetic field we estimated the syn- 
chrotron intensity at 23 GHz using the WMAP polarization 
data. This model is described in the next section. 



4. Synchrotron polarization at 23 GHz 

4.1. Synchrotron polarization 

Similarly to the synchrotron intensity (see Eq. 4) the stokes 
parameters Q and U of polarized synchrotron emission in- 
tegrated along the line of sight are 



Q{v) = fsf-s{v) / rie-B^^ ''^ COS 2(/) sin a rf^; 



(12) 



U{v) = fs^s{v) / ne-B^^'*'''''^^ sin2(/)sinad^; 

J z 



where 



cos = 



El ' 



Bl ' 



and 



sin a = v^l - Bl/B^. 



(13) 

(14) 
(15) 

(16) 



The polarization fraction fa is related to the cosmic ray 
energy distribution slope s: 



S+1 

s + 7/3' 



(17) 



For a typical value of cosmic ray spectrum s = 3 we obtain 
/. = 0.75. 

The polarized intensity is simply 



and the polarization angle 
7(;^) = 0.5tan-i (V{v) 



Qiy))' 



(18) 



(19) 



4.2. Polarization fraction at 23 GHz 

Combining equations 4, 12 and 13, the polarization fraction 
is given by the following simple equation: 



Siu) 



fs9 



(20) 



where fg is the intrinsic polarization fraction related to the 
cosmic rays energy spectrum (see Eq. 17) and g is a ge- 
ometrical reduction factor that reflects the depolarization 
due to the variation of the angle between the magnetic field 
and the line of sight directions. The g factor depends on the 
structure of the magnetic field and on the distribution of 
the cosmic rays in three-dimensions (i.e. the integrals in 
equations 4, 12 and 13). It takes values between and 1; it 
reaches a maximum when the magnetic field is parallel to 
the plane of the sky (B = B± ) . In practice it never reaches 
1 as the magnetic field direction fiuctuates on any line of 
sight, due to both the spiral structure of the magnetic field 
and turbulent motions. 

We can estimate the polarization fraction at 23 GHz 
using the WMAP polarization data (P23) and estimates of 
the synchrotron intensity (^23) given by models 1 to 3. The 
corresponding polarization fraction maps, smoothed at 5° 
are shown in Fig. 5. The polarization fraction for model 1 
is similar to the one given by Kogut et al. (2007) (their 
Fig. 4). Models 2 and 3 gives very similar polarization frac- 
tion maps, which is expected considering that both models 
have almost the same average f3s and considering the small 
dispersion of /3s in model 3. On the other hand the polar- 
ization fraction of model 1 is systematically lower than in 
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Fig. 5. Polarization fraction at 5° resolution. The polarization 
P is from the WMAP K band data and the intensity S is from 
models 1 to 3 (from top to bottom). 



other models because of the assumption of no anomalous 
emission which increases the synchrotron intensity (and not 
its polarization). Apart from the average level, the most 
significant difference between the polarization fraction of 
model 1 to 3 is the large scale structure. Models 2 and 3 
shows less polarization in the inner galaxy than towards 
the anti-center region. Model 1 has a systematically lower 
polarization value in the plane over all longitudes without 
a clear raise of the polarization fraction towards the anti- 
center. The contrast between low and mid latitudes is also 
stronger in model 1 than in models 2 and 3. 

The structures seen in the polarization fraction maps 
are due to real depolarization structure (due to the mag- 
netic field 3D structure - large scale and turbulent compo- 
nents) but also to our imperfect estimate of f3s ■ Nevertheless 
it is remarkable that the polarization fraction given by mod- 
els 2 and 3 are very similar even with important differences 
in the definition of (3s- In the following we look for a mag- 



netic field model that allows to reproduce these polarization 
fraction maps. 

4.3. Model of polarization at 23 GHz 

In model 4 we use the relation between intensity and po- 
larization to estimate the synchrotron intensity from the 
23 GHz polarization data which are assumed to be domi- 
nated by synchrotron polarization: 



5. 



23 



P23 

gfs ' 



The synchrotron spectral index is then given by: 



(3s = 0.248 log 



23 



gfsS. 



408 



(21) 



(22) 



This relies on a proper estimate of the reduction factor 
g. The construction of g depends on a model of the large 
scale galactic magnetic field and of its fluctuations due to 
turbulence. The latter will have a strong impact on the 
depolarization value, especially in the Galactic plane where 
lines of sight are longer. 

In order to estimate g, we use the fact that g ^ P/ S {fs 
is almost constant - see Eq. 17). To estimate the polariza- 
tion fraction P/S we use equations 4, 12, 13 and 18 which 
depends on the structure of the magnetic field B in three 
dimensions and on the electronic density rig . Following Page 
et al. (2007) we used a variation of the cosmic ray electron 
density with galacto-centric radius and height: 



= riQ exp(— r//ir) scch^(z/hz) 



(23) 



with hr — 5 kpc and = I kpc. Sun et al. (2007) adopted 
a similar description of the cosmic ray distribution. 

4.3.1. Galactic magnetic field at large scales 

Following what is found in the analysis of pulsar rotation 
measures (RM) (Han et al., 2006) we modeled the Galactic 
magnetic as a bi-symmetrical spiral (BSS): 



B{r,e,z) = Bo(r)cos(6l- Vln— )cos(x) 



(24) 



where ip = l/tan(p), p being the pitch angle and x(^) — 
Xotanh(z/l kpc) describing the vertical (z) component of 
the magnetic field (see Page et al., 2007). Using RMs at 
high Galactic latitudes, Han et al. (2006) estimated that 
the amplitude of this component of the magnetic field is 
almost 10 times smaller than the horizontal components. 
We found that xo has little impact on the modehng. 

Rotation measures are compatible with a slight varia- 
tion of the magnetic field strength Bq with Galactic radius 
R, modeled by Han et al. (2006) with a slowly rising ex- 
ponential towards the Galactic center. In fact this increase 
is not very clear and the data are also compatible with Bq 
constant with R. In this analysis we found that a constant 
Bq provides a better fit to the WMAP data. 

4.3.2. Turbulent part of the magnetic field 

The large scale spiral structure of the Galactic magnetic 
field is disrupted locally by turbulent motions of the gas. 
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Fig. 6. Polarization fraction as a function of depth along the 
line of sight. The large scale magnetic field is perpendicular 
to the line of sight, giving maximum polarization = 1 and 
fs = 0.75) with no turbulence. The different curves give the 
level of polarization we would observe for increasing level of the 
turbulent component of the magnetic field. 



The amplitude of the perturbation will depend on the tur- 
bulence energy density; on average data indicates that the 
energy in the turbulent part of the magnetic field is of 
the same order as the energy in the large scale field (see 
Beck (2001) and references therein). In addition Minter & 
Spangler (1996) showed that the power spectrum of the tur- 
bulent part of the magnetic field follows a power law with 
an exponent of -5/3 at scales smaller than 100 pc. At larger 
scales, the structure of the magnetic field is dominated by 
the spiral pattern. 

The impact of this turbulent part of the magnetic field 
on polarization is to lower the polarization fraction with 
respect to what would be observed with only the large scale 
magnetic field. This turbulent depolarization is stronger at 
low Galactic latitudes as lines of sight in the Galactic plane 
cross more fluctuations than at high latitudes. This effect 
is shown in Fig. 6 where we have simulated the effect of a 
turbulent magnetic field superposed on a large scale field 
perpendicular to the line of sight {B^ = By = 1 /iG and 
Bz = fJ.G). The x, y and z components of the turbulent 
field are independent ID Gaussian random fields along the 
line of sight (10 000 points, cell size of 0.1 pc) with a power 
spectrum following a -5/3 power law at scales lower than 
100 pc. For this exercise we assumed a constant cosmic ray 
density and s = 3 (i.e fs = 0.75). The polarization fraction 
is then computed by integrating numerically equations 4, 
12 and 13 for increasing values of the depth L. This was 
repeated for increasing values of Bturb from 0.01 to 10 fiG. 

As expected the polarization fraction decreases gener- 
ally with the depth of the line of sight and the amplitude 
of the turbulent motions. More and/or stronger turbulent 
fluctuations reduce the coherence of the field and the po- 
larization fraction. In fact the behavior observed here de- 
pends solely on the ratio of Bturb/ B and not on their re- 
spective absolute values. We noted that for moderate values 
of Bturb/ B, the polarization fraction reaches an asymptotic 
value. This values is reached at larger line of sight depths 
for greater Bturb/ B values. 
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Fig. 7. Model of the polarization fraction with (bottom) and 
without (top) turbulence. Parameters of the magnetic field 
model are p — —8.5°, Xo = 8°, tq = 11 kpc, Bq — 3 /iG and 
turbulence of aturt = 1.7 /iG. 



The map of g was computed by adding a turbulent vec- 
tor field to the BSS magnetic field model for each sky coor- 
dinate. Therefore every lines of sight are independent and 
no coherence on the plane of the sky due to turbulence is 
introduced. The turbulent part added here allows to set 
the proper average level of g but it also modifies the struc- 
ture of g as the depolarization due to turbulence depends 
on the length of the line of sight and the ratio of B±/B. 
Nevertheless most of the structure of g is set by the spiral 
structure of the field. 

An example of a gfs map, with and without turbulence 
added is shown in Fig. 7. Depolarization due to turbulence 
is increased in the galactic plane (where line of sight are 
longer) and in regions where B±/B is lower. The two bright 
spots above and under the Galactic center corresponds to 
regions where the Galactic magnetic field is almost perpen- 
dicular to the line of sight . 

4.3.3. Fitting the magnetic field parameters 

Unlike Page et al. (2007) we did not fit the parameters of 
the model on the polarization angle which is very sensitive 
to high latitudes structures like Loop 1. Instead the parame- 
ters are determined by fitting the polarization fraction map 
directly. Our model is not depedant on the absolute values 
of Bq and Bturb but only to their ratio. The parameters 
we are varying to fit the data are : p, tq, xo and Bturb- 
The fixed parameters are : jturb = —5/3, Imax = 100 pc, 
Bq = 3 /iG and Ue- 

The fitting of the parameters is done by minimizing 
the of the polarization fraction of our Galactic mag- 
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Fig. 8. Polarization fraction smoothed at 20°. Top: Model 
2: constant [3s = —3. Middle: Model 3: Ps estimated with a 
anomalous emission proportional to E[B — V). Bottom: Model 
4: Ps estimated from polarization for p = —8.5°, xo ~ 8°, 
ro — 11 kpc, Bq — 3 fiG and turbulence of aturt = 1-7 /iG. 



netic field model with respect to the polarization fraction 
deduced from models 2 and 3 and smoothed at 20° (see 
Fig. 8). Similar values are found for both model 2 and 3. 
The best fit model gives p = -8.5°, x = 8°, tq = 11.0 kpc, 
Bturb/Bo = 0.57. The polarization fraction of our best fit 
model is shown in Fig. 8. The correlation factor are 0.65 
and 0.72 for model 2 and 3 respectively. 

The parameters found for the large scale field are in very 
good agreement with what is obtained from the analysis of 
pulsar rotation measures and from a previous analysis of 
the WMAP data by Jansson et al. (2007). In particular 
the pitch angle value we found (p = —8.5°) is remarkably 
close to the values obtained using pulsar rotation measures 
(Han & Qiao, 1994; Han et al., 1999; Indrani & Deshpande, 
1999; Beck, 2001; Han et al., 2006) or polarized starhght 
absorption measurements (Heiles et al., 1996) where values 
are in found in the range —7 < p < —11°. On the other 




Fig. 9. Top: Polarized intensity (in mK) at 23 GHz, smoothed 
at 5°. Bottom: Mask showing regions of signal-to-noise ratio 
greater than 3. 



hand our result is not compatible with the result of Page 
et al. (2007) who obtained p = -35°. 

The values of Bturb / -So found for model 2 and 3 is also 
remarkably close to one obtained from analysis of the syn- 
chrotron emission (Phillipps et al., 1981; Beck, 2001) which 
gives Bturb/Bo ~ 0.66, which is implies energy equipar- 
tition between the turbulent and large scale components. 
Our result is also in agreement with the work of Jones 
et al. (1992) who obtained Bturb/ Bq ~= 0.6 using po- 
larized extinction measurements at 2.2 ^m. On the other 
hand, and similarly to Sun et al. (2007) who did not con- 
sider the presence of anomalous emission at 23 GHz, model 
1 requires Bturb/Bo ^ 1.5 to fit the data, which is signifi- 
cantly stronger than the value obtained by Phillipps et al. 
(1981); Beck (2001). This is due to the fact that anoma- 
lous emission contributes a significant fraction of the in- 
tensity at 23 GHz and that it is not polarized. If (unpo- 
larized) anomalous emission is not taken into account at 
23 GHz, a stronger depolarization (and therefore more tur- 
bulence) is needed to account for the low polarization frac- 
tion. The incapacity of model 1 to reproduce the right level 
of Bturb / Bq IS & strong evidence for the presence of an un- 
polarized anomalous emission at 23 GHz. Therefore in the 
following we will not consider model 1 anymore. 

4.4. Model of synchrotron emission 

We build the synchrotron spectral index using 

. _ log(M23/^408) 

~ log(23/0.408) ^ > 



Miville-Deschenes, M.-A. et al.: Anomalous and synchrotron emissions using WMAP polarization data 



9 



where M23 is the model of the intensity at 23 GHz obtained 
from polarization: 

M23 = ^. (26) 

gjs 

P23 is the polarization map from WMAP smoothed at 5° 
(shown in Fig. 9), fs is given by Eq. 17 and the depolar- 
ization fraction g is given by our model, smoothed at 20° 
to remove any small scale fluctuations introduced by the 
turbulent part of the magnetic field. 

A reliable estimate of M23 can only be made for pixels 
where the signal-to-noise ratio of the polarization data is 
significant. We defined a mask shown in Fig. 9 which gives 
regions with polarization signal-to-noise ratio large enough 
to estimate the synchrotron intensity. Inside this mask (3s 
was estimated using Eq. 25. Outside the polarization mask, 
E(B-V) is mostly lower than 1 which corresponds to regions 
of the sky where the anomalous emission correlates well 
with extinction. Therefore we used the estimate of Ps given 
by model 3 in regions outside the polarization mask. The 
resulting /3s map at a resolution of 5° is shown in the bottom 
panel of Fig. 2. Even if Ps was estimated with two different 
models, we note that there are no striking discontinuities 
at the boundaries of the mask. Compared to models 1 and 

3 the galactic plane does not stand out in this spectral 
index map. In that respect the spectral index map deduced 
from the polarization data is much more similar to the ones 
determined between 408, 1420 and 2326 MHz by Giardino 
et al. (2002) and by Platania et al. (2003). 

The histogram of j3s is shown in Fig. 3 (dashed line). 
It has an average value of -3.00 and a standard deviation 
of 0.06. The minimum dispersion of Ps is thus obtained for 
model 4. 

5. Anomalous emission in the WMAP data 

The predicted synchrotron intensity at 23 GHz for mod- 
els 2, 3 and 4 is shown in Fig. 10. For model 3 (Fig. 10- 
middle), there is more synchrotron emission in the plane 
than for the two other models. This is due to the spectral 
index used which is significantly shallower in the plane (see 
Fig 2-middle). The synchrotron emission for models 2 and 

4 are very similar in overall intensity but there are signif- 
icant spatial variations, especially in the Galactic plane, 
partly due to the fact that model 4 does not include spatial 
variations of g at scales smaller than 20°. 

The residual emission, when free-free and synchrotron 
are removed from the 23 GHz data, are shown in Fig. 11 for 
models 2, 3 and 4. We attribute this residual to the anoma- 
lous emission. The anomalous emission obtained for models 
2 (constant (3s — —3) and model 4 (/3s obtained from po- 
larization) are very similar. The difference between these 
two models is well fitted by a Gaussian of average 0.002 mk 
and a = 0.015 mK. For both models, the correlation fac- 
tor of the anomalous emission with extinction (E(B-V)) is 
0.94. This analysis shows that the synchrotron emission at 
23 GHz can be rather well estimated, to first approxima- 
tion, with an extrapolation of the 408 MHz with a constant 
spectral index Ps — —3. In the following we describe the 
properties of the anomalous emission for model 4. 

The anomalous emission is dominating the signal from 
the Galactic plane to mid-latitude regions, with a mean in- 
tensity up to ten times higher than synchrotron (see Fig. 




Fig. 10. Synchrotron emission (in logio(mK)) at 23 GHz for 
models 2 (top), 3 (middle) and 4 (bottom). The synchrotron 
map for model 2 (top) is a scaling of the 408 MHz map with a 
constant (is ~ —3. 



12). On the other hand, we note that the average level of 
synchrotron emission at high latitude is about three times 
that of the anomalous emission, likely due to the large scale 
height of cosmic ray electrons. The Gould belt's cloud have 
strong anomalous emission and known high-latitude cir- 
rus regions like the North Celestial Loop are also seen in 
anomalous emission. 

The anomalous emission was also extracted at 33, 41 
and 61 GHz using the synchrotron spectral index map of 
model 4. The emissivity per hydrogen atom at each fre- 
quency of the anomalous emission was estimated by corre- 
lating it with Nh (cm-2) = 5. 8 x 10'^^ E{B-V) (see Fig 13) 
for Nh > 10^^ cm~^. The spectrum obtained using model 4 
is shown in Fig. 13; the anomalous emission spectrum is the 
same within 1% for models 2 and 4. In Fig. 13 the anoma- 
lous emission spectrum is compared to the spinning dust 
model of Draine & Lazarian (1998b) for the Warm Neutral 
Medium ( WNM) and with the model of Ysard & Verstraete 
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Fig. 12. Ratio of anomalous emission to synchrotron at 23 
GHz (model 4) as a function of column density. In the plane, 
the anomalous is on average several times higher than the syn- 
chrotron emission. In the most diffuse regions at high Galactic 
latitude the synchrotron emission overcomes the anomalous 
emission in total intensity. The error bars represent the stan- 
dard deviation of the anomalous to synchrotron ratio in each 
bin of column density. 




Fig. 11. Estimate of the anomalous emission (in logio(mK)) at 
23 GHz for models 2 (top), 3 (middle) and 4 (bottom). 




(2008) (for PAH sizes ranging from 50 to 96 carbon atoms). 
Both models reproduce well the observed spectrum. A de- 
tailed analysis of the anomalous emission spectrum is left 
to a further study. 

6. Conclusion 

In this paper we presented a new analysis of the Galactic 
emission at 23 GHz as observed by WMAP. The goal of this 
study was to combine the total intensity and polarization 
data from WMAP to separate the synchrotron emission and 
the anomalous emission at this frequency. Our analysis is 
based on the hypothesis that 1) the 23 GHz intensity data 
of WMAP is dominated by CMB, synchrotron, free-free 
and anomalous emission, 2) the 23 GHz polarization data 
is dominated by synchrotron only 3) the average spectral 
index of synchrotron between 408 MHz and 23 GHz is close 
to = —3, and 4) the Galactic magnetic field has a BSS 
large scale structure with a turbulent part that follows a 
-5/3 power spectrum at scales smaller than 100 pc. 



Fig. 13. Spectrum of the anomalous emissivity obtained by 
correlating the anomalous from model 4 with A^'^f = 8.0 x 
10^^£(B - V). The correlation was estimated for Nh > 
10^^ cm~'^. The uncertainty is the standard deviation of the 
emissivity observed for the selected pixels. The emissivity at 
61 GHz is compatible with zero. The curves are the emissivity 
spectra, normalized at 23 GHz, given by Draine & Lazarian 
(1998b) for the WNM (dashed) and by Ysard & Verstraete 
(2008) for PAHs ranging from 50 to 96 carbon atoms (dotted). 



We used the WMAP 23 GHz polarization data to con- 
strain the parameters of the magnetic field. The data are 
found compatible with a pitch angle of —8.5° and a ratio 
of the turbulent to regular part of the magnetic field of 
Bturb/Bo = 0.57. The parameters of our Galactic magnetic 
field model are in good agreement with what is found with 
pulsar rotation measures, radio synchrotron emission and 
polarized extinction in the near-infrared. 

Based on our Galactic magnetic field model and the 
WMAP polarization data, we produced a map of the syn- 
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chrotron spectral index in total intensity, between 408 MHz 
and 23 GHz. The map of fig we found does not show a strong 
Galactic plane feature, in accordance with what was found 
in the radio by Giardino ct al. (2002); Platania ct al. (2003). 
Extrapolating the Haslam 408 MHz with this new estimate 
of the synchrotron spectral index, we could separate the 
contribution of the synchrotron and anomalous emission at 
23, 33, 41 and 61 GHz. 

The all-sky anomalous emission at 23 GHz is found to 
be correlated to E(B-V) (or to Ay) with a correlation fac- 
tor of 94%. It is well detected in the Galactic plane, in 
Gould Belt's clouds and in the main high-latitude cirrus 
clouds. In the Galactic plane the anomalous emission is up 
to ten times strongc^r than the synchrotron but it is only one 
third of the synchrotron in the most diffuse regions of the 
sky where the synchrotron, while very smooth, has an av- 
erage level stronger than the anomalous emission, possibly 
due to the larger scale height of cosmic ray electrons com- 
pare to interstellar dust. Under the hypothesis that the syn- 
chrotron spectral index does not vary with frequency over 
the WMAP range, we found that the anomalous emission 
has a spectrum from 23 to 61 GHz in accordance with the 
model of spinning dust from Draine & Lazarian (1998b). A 
deeper analysis of the origin of the anomalous emission and 
of its consequences regarding the properties of interstellar 
dust will be the subject of a forthcoming paper. 

The presence in the 20-200 GHz range of several diffuse 
Galactic components with similar intensities and some spa- 
tial correlation makes the extraction of the CMB a complex 
task. In order to achieve sufficient accuracy on the cos- 
mological signal the component separation needs to take 
advantage of the knowledge on the properties of diffuse 
Galactic emission. In this paper we have shown that com- 
bining temperature and polarization data helps to improve 
our model of Galactic emissions. The higher signal-to-noise 
ratio, the higher angular resolution, the larger frequency 
coverage and the increase number of frequency bands of 
the Planck data will also certainly help to make progress in 
this area. 

Some of the results in this paper have been obtained 

using the HEALPix package (Gorski et al., 2005). This work 
was partly supported by the Canadian Space Agency. 
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